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Backlund Transformations of Soliton Systems from 
Symmetry Constraints 

Wen-Xiu Ma and Xianguo Geng 



Abstract. Binary symmetry constraints are applied to constructing Backlund 
transformations of soliton systems, both continuous and discrete. Construc- 
tion of solutions to soliton systems is split into finding solutions to lower- 
dimensional Liouville integrable systems, which also paves a way for separa- 
tion of variables and exhibits integrability by quadratures for soliton systems. 
Kyj Illustrative examples are provided for the KdV equation, the AKNS system of 

_^ , nonlinear Schrodinger equations, the Toda lattice, and the Langmuir lattice. 



1. Introduction 

Symmetry constraints w, B play an important role in showing integrability 
by quadratures for soliton systems, both continuous and discrete. They help to 
generate finite-dimensional integrable systems pL pL and integrable symplec- 
tic mappings pi , and further provide a way of constructing finite-gap solutions to 
f^ , soliton systems by means of Riemann-theta functions H]. Based on Lax pairs, sym- 

metries themselves can also be applied to constructing Backlund transformations 
of soliton systems m . However, there is little work showing the importance of sym- 
metry constraints in the study of Backlund transformations. In this paper we focus 
on the construction of Backlund transformations by using symmetry constraints, 
C ' and show in certain cases that symmetry constraints can break up soliton systems 

^ , into lower-dimensional Liouville integrable systems. 

Let us recall some fundamental concepts. A system of continuous equations 
ut — K{u, Ux, ■ ■ ■) is said to have a continuous Lax pair 



X 

. .^ , (1-1) (/., = C/(u, A)./., (^t = y(M, u„ • • • ; A)(^, 



if it is equivalent to the compatibility condition Ut — Vx + [U,V] = of (1.1) 
under the isospectral condition Aj = 0, and a system of discrete equations ut = 
K{u, E~^u, Eu, ■ ■ ■) (where E is the shift operator, i.e., {Eu){n) ~ u(n + 1)) is said 
to have a discrete Lax pair 

(1.2) Ecj) = U(u, A)0, (t>t = V{u, E-^u, Eu,--- ; X)(j), 
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if it is equivalent to the compatibility condition Ut — {EV)U+UV = of (1.2) under 
the isospectral condition Af — 0. The corresponding adjoint Lax pairs presenting 
the same compatibility conditions read as 

(1.3) ^. = ~U^{u, X)4>, ^t = -V^{u, u,,---; A)Vj; 

(1.4) ^- V = {E-^U'^{u, A))V', ^t = -V^iu, E'^u, Eu,--- ] A)V'; 

where (•)"^ denotes matrix transpose. Adjoint Lax pairs can help us determine the 
variational derivative of the spectral parameter with respect to the potential u (see, 
for example, |g, |^ for the continuous case). 

A soliton hierarchy of continuous or discrete systems 

(1.5) Ut,, = i^„ - $"i^o - JGn = J— ^, n > 0, 

ou 

can be generated through the isospectral (At„ = 0) compatability conditions of 
continuous Lax pairs 

(1.6) (l>,^U{u,X)<l>, </.«„= y(")(u,7.,,-.-;A)0, 
or discrete Lax pairs 

(1.7) E(l) = U{u, A)</>, (/)t„ = y("^ (u, E-'^u, Eu,--- \ A)(/>, 

where F^") are Laurent polynomials in A, <& is a hereditary recursion operator to 
map symmetries into symmetries, and J is a Hamiltonian operator to map conserved 
covariants to symmetries. 

In this paper, we would like to show that symmetry constraints can be applied 
to constructing Backlund transformations of soliton systems from Lax pairs. The 
resulting Backlund transformations separate each soliton system (in a hierarchy) 
into two lower-dimensional integrable systems. Thus, symmetry constraints are 
shown to be very useful in exposing integrability by quadratures for soliton systems. 
Illustrative examples will be given in both continuous and discrete cases. 

2. Symmetry constraints 

Let us consider the space parts and the time parts of Lax pairs and adjoint Lax 
pairs: 

r 0, = C/(m, A)0, r 0t„ = !/(")(«, u,, • • ■ ; A)0, 

1 ^. = -f/^(«, A)V; 1 ^t„ = -y(»)^(u, «„•••; A)^; 

or 

f £;</) = C/(u, A)0, J 0t„ =y(")(u, £;-!«, £;u,--- ;A)0, 

\ E-^ij = {E-^U^{u, X))tP; \ TPt„ = -y(")^(u, E-'^u, Eu,--- ] A)V'. 

By using the space parts, we can work out 

(2.1) T- = «V' ^., 4^. or — = /3(^V 



5u du ' 5u du 

5u 



where a and (3 are two constants. Note that the Lie homomorphism J 4- transforms 
conserved functionals to symmetries. Therefore, j|^ is a symmetry of each system 
Ut„ — Kn, since A is a conserved functional, i.e., At„ — {X{u))t,, = when ut,, — Kn- 
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Now for all ttiq > 0, we can make symmetry constraints 

(2.2) K„,^ = -J— = J^^ \ V , or K„,, = -J— = J{E4,^) ^0. 

If wc take distinct eigenvalues Ai, • • • , Ajv, and suppose that 

(2.3) (t)^^^ = C/(w, \j)cl)^^\ V^i^' = -[/^(u, \j)tp^'\ 
or 

(2.4) E(l)^i^ = U{u,\j)(j,^^\ E-^t/j^^^ == iE-^U^{u,Xj))tlj^^\ 
where 

0(^') = (01„ • • • ,<j>sjf, ^(^'^ = (V'l,, • • • ,iJsjf, 

we can make more systematical symmetry constraints 

namely, 

(2.5) X™, = f: J^W^^^^I^^W, or X„„ = f: J(i^^O-)^)^^^^<^0-), 



where aj and /3j are the constants defined as in ( 2A ) . These symmetry constraints 
suggest 

(2.6) G,„„ = f: V«^^^^|^</>'^"\ or G^o = ^(iJV.'^")^)^^^!^^^^'. 

Among those constraints between the potential, u, and the eigenfunctions and ad- 
joint eigenfunctions, (j)^^' and ip , the Bargmann constraint will be applied to con- 
structing Backlund transformations between soliton systems and lower-dimensional 
integrable systems. 

3. Backlund transformations 



Let us take the Bargmann constraint, i.e., the constraint (2.6) with Gmo = 
Gmoiu) not involving any d^u {d — ■^), i > 0, or any E^u, i =^ 0. Note that the 



discrete constraint defined by (2.6) can be rewritten as 



N , r.. N 



(3.1) G„M = %^ =T.T-^ = Y.^^^^-U-\u,X,fJ^^^^^\ 

ou ^-^ Pi ou ^-^ ou 

and therefore, in each of continuous and discrete cases, the Bargmann constraint 
defined by ( p^ ) is an algebraic equation on m, 4>^-^\ and ■0'^.'). Assume that solving 
the corresponding algebraic equation for u gives rise to an explicit expression of u: 

(3.2) . = /(0W,0(2),...,0W;V^W,^(^),...,^W). 

Substituting this expression of u into Lax pairs and adjoint Lax pairs leads to two 
systems, called binary constrained Lax pairs. Binary constrained continuous Lax 



pairs read as 
(3.3) 

(3.4) 
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# = [/(/, Aj)0(^'\ 1<J<^, 
4^"U-t/^(/,Aj)^«, 1<3<N; 



the first of which is a system of ordinary differential equations, but the second of 
which is a system of partial differential equations since V^"^' contains some deriva- 
tives of u with respect to x. Binary constrained discrete Lax pairs read as 

r E<j>^i) - [/(/, A,)(/)(J") , 1 < j < TV, 

I E-^i^(i) = {E-'U^if,X,))ij(j\ l<j<N- 

, , \ 4^1 - ^("H/, E-^!. Ef, ■■■; A,)0O"), 1 < J < TV, 

(3.6) { 

[ ^i^' = -y(")^(/,i?-i/,i?/,---;A,)^(^'\ i<j<A^; 

the first of which is a system of difference equations, but the second of which 
is a system of difference-differential equations since y'"-' contains some of E^u, 
i ^ 0. However, the second systems can be tranfornied into systems of ordinary 
differential equations by using the first systems. Furthermore, it can be shown by 
r-matrix formulation that all binary constrained Lax pairs, both continuous and 



discrete, are integrable in the Liouville sense |1C, 11 



Therefore, (3.2) provides Backlund transformations between soliton systems 
and integrable binary constrained Lax pairs, and construction of solutions u = 
f{(l)^^\ip^^^) to soliton systems is split into finding solutions t/)*^-'^ and i/'^"''' to two 
lower-dimensional integrable systems. 

4. Examples of continuous systems 

Example 4.1. Let us consider the KdV Equation 

1 , 3 JHi , ^ ~, fA ,33 



(4.1) uti = -u^xxx + 7:""^ = J^—^ J ^ d, Hi = / {-uu^j: + -—u ) dx, 

4 z ou J a 12 

which can be written as Ut — Vx + [U, V] — with 
(4.2) 




2""' 4 
Take the Bargmann symmetry constraint 

N 

du 



K,^„ = 2dx J2 V,0)^ ^^("'^^V ^^\ where K^„ ^ K, = u. 



which implies the following equation 

^ / \ / 01, \ J^ 
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Integrating this equation with respect to x yields a Backhmd transformation 

N 
(4.3) U = f{(|)^f,i^^,) = 2Y,<l^l3i>2j - 2 < Pi,Q2 >, 



J = l 



,JV 



where < ■, ■ > denotes the standard inner product of R , and Pi and Qi are defined 

by 



Pr 



IJil,(Pi2, ■ ■ ■ , <PiN) , ^i 



i'iptl,A: 



, V'iw)^, i = 1,2. 



A general Backlund transformation u ~ 2 < Pi,Q2 > +c with an arbitrary constant 
c can also be resulted from the above symmetry constraint, but it will not generate 
essentially new integrable systems from Lax pair of the KdV equation and so we 
omit to discuss it. 

Keeping (4.3) in mind, the corresponding constrained Lax pairs (3.3) and (3.4), 
where two matrices U and V'^"^ — V are defined by (4.2), can simultaneously be 
rewritten as 



(4.4) 
(4.5) 



P..r = -^^, 0„ = -T^, H- = -Fs, I = 1,2, 



Pit.^ — 



dHI 



Qiti — 



dP^ 
dHI 



dQ,' ^"^ dP, 
where the functions Fm are given by 
Fo = 0, Fi = 1, F2 = 0, 

F3^<P2,Ql > -<Pl,Q2 >^ + <APi,Q2 >, 



H'i = -Fi, z = l,2, 



F,n=Y.^ 



aiam-i-A + biC,n-i-i) + C„i-3 + bm-2 - x/&m-3, TO > 4, 



a, =< A'Pi,Qi >-< A'P2,Q2 >, k =< A'Pi, O2 >, c, =< A'Fa, Qi >, 

through F = detV\u=f = E™=o^™^"'"' ^ = E»=o ^^^"' satisfying 14 = [U,V]. 
Note that we always accept 



A^ diag(Ai,--- ,A 



N 



These two systems ( |4.4| ) and (4^) are Liouville integrable |1C], since they have 
involutive integrals of motion Fm, m. > 0, and 



(4.6) 



Fj = (jilji^lj + 4>2j'tp2j, ^ < j < N, 



among which F3, F4, • • • , Fm+2 and Fi, • • • , F/v are functionally independent. All 
functions 



X ti 



(4.7) u{x,h) = 2 < Pi{x,ti),Q2{x,ti) >=< gH9H,Pw.9H9H,Q20 > 



where g|f and g^ are the Hamiltonian flows associated with (4.4) and (4.5) respec- 
tively, and Pio and Q20 are two arbitrary constant vectors, will determine solutions 
to the KdV equation (4.1). 



(4.8) 



Example 4.2. Let us now consider the coupled nonlinear Schrodinger system: 
\ r 



Ut2 



t2 
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with the Hamiltonian operator J and the Hamiltonian functional H2 
( 



(4.9) J = I 

It has a Lax pair 



, -f^2 = T^ / {qr^^ - q^r^ + q^^r - iq^r^) dx. 



(4.10) 



U = 



-A 



A 



, v = 



-A^ + igrr q\ - \q^ 
rX + irj; 



A^ - ^qr 



The Bargmann symmetry constraint (2.5) reads as 



K^^ ^ jf2i^U)TdU(^^U)^ ^here if„„ = i^o = J ( ^ ) 
j=i \ / 



This imphes the foUowing equation 



J 



/ / 1 




which equivalently engenders a Backlund transformation 

N 



N 



(V'lj,V'2j) 
(^lj,V'2i) 



01J 
02j 

02j 



N 



(/'2jVlj 
(/'ljl/'2j 



(4.11) 



/(^.;^.)-E(^::J: 



j=i 



<PuQ2> 

<P2,Ql > 



Keeping ( 4.11 ) in mind, the corresponding constrained Lax pairs (3.3) and (3.4), 
where two matrices U and F*^"^ — V are defined by (4.10), can simultaneously be 
transformed into the following 



(4.12) P„ 

(4.13) Pu, 



dm 



, Q^. 



?^, H-^F2--Fl z = l,2, 
dp ^ 4 1 ' 



dm 



dQ.r "" dP, 
where the functions Fm are given by 

Fi=<P2,Q2>-<Pl,Ql>, 

in -2 



m = F. 



-F1F2 + —F?. i = 1,2, 



F,n = 51 [i(< ^'-Pl'Ql > - < ^'-P2,Q2 >)(< A'"-*-2Pl,Ql > 



< A™-*-2p2, Q2 >)+ < A'Pi, Q2 X A"'-'-^P2,Qi > 

< A'"-1P2, Q2 > - < A™"'Pi, Qi >, m > 2. 



These two systems are completely integrable in the Liouville sense [10|, since they 
have involutive integrals of motion Fm, to > 0, defined above, and Fj, 1 < j < N, 
defined by (16), among which Pi, F2, ■ ' ■ , P/v and Pi, ■ ' ■ j Pv are functionally inde- 
pendent. The Backlund transformation (111) determines solutions to the coupled 
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nonlinear Schrodinger system (4.8): 

q{x,t2) =< Pl{x,t2),Q2{x,t2) >= y^^<j)ij{x,t2)lp2j{x,t2), 



(4.14) 



N 



3 = i 

N 



"(X,t2) =< P2{x,t2),Ql{x,t2) >= y^^(j>2j{x,t2)^plj{x,t2), 



if (I)ij{x,t2) and TJJij{x ,t2) simultaneously solve two integrable finite-dimensional 
Hamiltonian systems ( ff.l2 ) and ( 4.13 ). 

5. Examples of discrete systems 

Example 5.1. Let us consider the Toda lattice Jl^ : 

(5.1) at{n,t) — a{n,t){b{n + l,t) — b{n,t)), bt{n,t) = a{n,t) ~ a{n ~ l,t), 
which associates with the discrete spectral problem 

(5.2) Ecf>^U{u,X)q^, U{u,X)^l ^^ a"& )' "^^ ( ^2 )' 

where u = (a, 6)"^ and A is a spectral parameter. In order to derive a hierarchy of 
lattice equations associated with ( |5.2| ), we first solve the stationary discrete zero- 
curvature equation: 

(5.3) {EV)U-UV^O, V^{V,,)2y,2, 
by assuming 

Vn =aB + {b- X)C, Vu = E-^aC, 1^2i = -C, ^22 = E-^aB, 
where 



The discrete zero-curvature equation (5.3) requires 

(5.4) JG_i = 0, MGn-i = JGn, n>0. 

where G„ = {Bn,Cn) , and J, M are two skew-symmetric operators: 



(5.5) 



J = 



aA \ ( a(A-A*)a aAb \ 

-A*a / V -bA*a aA-A*aj' 

We choose G_i = (0, 1)"^, and assume that all terms of G„, n > 0, do not belong 
to kerJ = span{G_i, G_2} where G_2 = {ci~^,0) , when we uniquely determine 
all G„, n > 0. For instance, the second member has to be Go = (1,6)-'". This 
requirement also means that we just choose the key vector fields to form systems 
of lattice equations. 

Let Ai, • ■ • , Xn be distinct eigenvalues. Then we have 

(5.6) iE^ij,E<f>2j) = i<j>i,,h3)Uiu, X,f, 1<J<N; 

(5.7) iE4'i,,Ei,2j) - {^ij,^P2j)Uiu, Xj)-\ l<j<N. 
It is easy to see that 



(5.8) 



M. 



J_X^ 
5u 



= X,J 



SX^ 
Su 
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where 6Xj/Su is determined by 



(5.9) 



6u 



sx 



-d \ 13. ( (Aj - b)4>2jiJlj + 02iV'2j 



5Xj 
Sb 



a<j)2ji'i] 



, Pj = const. 



Now the Bargmann constraint Go = "^j^i Pj ^6\j/Su leads to a Backhmd trans- 
formation 

(5.10) a=<AP2,Pl> + <P2,Q2>-<P2,Ql>^ b=<P2,Ql>, 

where A = diag(Ai,--- ,Xn), Pi = (0ii,--- ,4'inV, Qi = (V'ii,"- ,'4'in)'^ , and 
< •, • > is the standard inner product of R^, as defined before. 

Substituting (5.1C) into ( [3.6| ) and (5.7) yields a discrete Bargmann system 

r EPi = (< AP2,Ql > + <P2,Q2>-< Pi, Qi >') P2, 
£^^^2 = -Pi- < P2, Ql> P2+ AP2, 



(5.11) 



EQi 



Q2- < P2,Qi > Qi + AQi 



< AP2,Qi > + < P2, Q2 > - < P2, Qi >2 ' 
L EQ2 = -Qi, 
which determines a symplectic mapping H: 
(5.12) {EP,,EP2,EQ,,EQ2)=H{P,,P2,Qi,Q2), 

since we have by a direct calculation 

AT AT 

J2 d{E(l)^^'>) A d{Etlj^^'>) = ^ d0(j') A #(^'\ 

The generating function J^a — det y|,j=/: 

^A = -QxiAPuQi) - Qx{Pi,Q'2)+ < PuQi > Qx{P2,Qi) 
Qx{Pi,Qi) QxiPuQi) 

Qx{P2,Qi) Qx{P2,Q2) 



= ^ F„A-"-\ 



m>0 



where 



AT 



Qa(C,^) = E 



J = l 



A~A, 



Y, <A™e,r,>A-'"-\ 



m>0 



generates a hierarchy of invariants of ( ^.11| ) : 

Fo=<APi,Qi >-<Pi,Q2> + <Pi,Qi><P2,Qi >, 
F™ = - < ^"+iPi, Qi > - < A'"Pi, Q2 > + < Pi, Qi >< A'"P2, Qi > 

< A'-iPi, Qi > < A™-'P2, Qi > 

< A'-iPi, Q2 > < A™-*P2, Q2 > 
A direct computation can show the involutivity 

{P„,P,} = 0, m,l>0, 

where the variants Fj = ^iji/'ij + 4>2j'^2j, 1 < j < ^, defined as before. Now we 
can easily see that the symplectic mapping (5.11) is Liouville integrable |11]. 



■E 



m > 1. 
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Introduce an initial-value problem 

(5.13) P't^g^, Qu = ~-g^, {P^it),Q^mt=0 = {P^0,Q^0): ^=1,2, 

where -P^o and Qio, i = 1,2, are arbitrary constant vectors. Let {Pi{t),Qi{t)), i = 
1,2, be a solution to the initial- value problems (5.13), and further define 

(5.14) (Pi(n,i),P2(n,t),Qi(n,i),Q2(n,t)) = i/"(Pi(i),P2(i),Qi(i),Q2(t)). 

Then a{n,t) and b{n,t) determined by the Backlund transformation (p3Q) solves 
the Toda lattice (O). 



Example 5.2. Let us now consider the Langmuir lattice |13|: 

(5.15) o,t{n, t) — a{n, t){a{n -|- 1, t) — a{n — 1, t)), 

which associates with a reduction of the discrete spectral problem ( |5.2| ): 

/ a \ / </>! \ 

(5.16) E(j)^U{a,\)(l), U{a,X)=\ \ cf) = \ 

Assume that its discrete zero-curvature equation 

(5.17) {EV)U-UV^{), F=(F,,)2x2 
has a solution 

Fii =aB- X^C, ^12 = XE-'^aC, V21 = -AC, ^22 = E-'^aB, 

where 

B = Y,B^X-^\ C = ^C.A-^\ 

i>0 i>0 



Then upon choosing Bq = Cq = 1, we can easily find that (5.17) is equivalent to 

(5.18) MBn-i = JB„, Cn = (A*a - aA)-^A*aBn, n > 1, 
where two skew-symmetric operators J and M read as 

(5.19) M^a{A-A*)a, J ^ aA{A*a - aA)-^A*a. 
Let Ai, • • ■ , Ajv be distinct eigenvalues, then we have 

r (^01, , £;02,) = (01, , 02,)(7(a, A,)^, 1 < J < ^, 
\ {Ei^i,,E^2j) = (^ij,iJ2j)U{a,X,)-\ 1<3<N, 
and 

(5.21) M—^ = A^ J— i-, where — ^ = ^(Aj(?!)2jVij + <t>2j'il}2j) , Pj = const. 

da ■'da da a 

Now similarly, the Bargmann constraint Go = Bq = X^^i 0j^ SXj/Sa leads to a 
Backlund transformation 

(5.22) a=<AP2,Qi> + <P2,Q2>- 



Substituting (5.22) into (p^2G), we obtain another discrete Bargmann system 
EPi = (< AP2,Qi > + < P2, Q2 >) P2, EP2 = -Pi + AP2, 

^^^ - <AP2,Qi> + <P2,Q2>' ""^^ - '^'' 
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which determines a symplectic mapping H: 

(5.24) iEPi,EP2,EQi,EQ2) = H{PuP2,Qi,Q2)- 



This symplectic mapping is Liouville integrable [11|, since we have involutive in- 



variants: Fj = (pijipij + 4'2j'4'2j, 1 < i < -^1 and Fm, m > 0, defined by 
Fo = - <A^Pi,Qi > -<APi,Q2 > + <Pi,Qi> i<AP2,Qi > + < P2,Q2 >), 

F™ = - < A2™+2Pi, Qi > - < A2"+lPi, Q2 > + < Pi, Ol >< ^2™ + lp2, Qi > 

< A'^'-'^Pi,Qi > < A^"'-^'+^P2,Qi > 

< A^'-^Pi,Q2 > < ^2™-2»+2p2, Q2 > 

Let {Pi{t), Qi{t)), i = 1, 2, be a sohition to an initial-value problem: 



< P2,Q2 >< A^"'Pi,Qi > +J2 



BF BF 

(5.25) ^«* = ^' ^"""af^' (^^W'Q^W)l*=o = (P^o,Q^o), i = l,2, 

where P^o and Qio^ i — 1,2, are arbitrary, and similarly define 

(5.26) {Pi{n,t), P2{n,t),Qi{n,t),Q2{n,t)) ^ H"iPi{t), P2{t),Qi{t),Q2{t)). 



Then a{n,t) determined by the Backlund transformation (5.22): 

a{n,t) =< AP2{n,t),Qi{n,t) > + < P2{n,t),Q2in,t) > 
provides a solution to the Langmuir lattice (^.15). 



6. Concluding remarks 

It has been shown that solving symmetry constraints for u can give rise to 
Backlund transformations between soliton systems and lower-dimensional Liouville 
integrable systems, which supplements the study of binary nonlinearization of Lax 
pairs 0, pL E4L 15 1 . Construction of solutions to soliton systems is split into 



finding solutions to the space and time parts of integrable constrained Lax pairs, 
which gives a way to separate variables for soliton systems and exhibits integrability 
by quadratures for soliton systems. Upon solving the Riemann-Jacobi inversion 
problems for constrained Lax pairs, the resulting Backlund transformations can 
generate finite-gap solutions to soliton systems in terms o f Rie mann-theta functions. 
We remark that all symmetry constraints defined by ( p.6[ ) can put Lax pairs into 
integrable symplectic mappings and/or integrable finite-dimensional Hamiltonian 
systems. The corresponding constrained Lax pairs may have some specific prop- 
erties, e.g., bi-Hamiltonian and quasi-bi-Hamiltonian structures. Therefore, sym- 
metry constraints are very powerful in constructing lower-dimensional integrable 
systems from Lax pairs of soliton systems. Nevertheless, there exist symmetry con- 
straints which do not force Lax pairs into integrable systems |16[| , and the problem 
of integrability has not been solved for the time parts of the original constrained 
Lax pairs, i.e., systems of partial differential equations 

,itf = ^^"^ (/,/.,•••; A,)0(^) , 1 < J < A^, 

^U) ^ _^(n)T(j^ f^^.... x^)^U)^ 1<J<N; 
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and systems of difFerence-difFerential equations 

, , f '/'if = ^^"H/, Ef, E-^f, • • • ; \,)<J>(J\ 1<J<N, 

(6.2) < 

[ 4'J - -V(")^(/, Ef, E-^f, ■■■; A,)^(^), 1 < J < iV. 

We are curious to known whether they are good candidates for integrable systems. 
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